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Abstract 

Explicit analytical expressions for Helmholtz free energy, chemical poten- 
tial, entropy and pressure of the multi-component dimerizing Yukawa hard- 
sphere fluid are presented. These expressions are written in terms of the 
Blum's scaling parameter F, which follows from the solution of the associative 
mean spherical approximation (AMSA) for the model with factorized Yukawa 
coefficients. In this case solution of the AMSA reduces to the solution of only 
one nonlinear algebraic equation for F. This feature enables the theory to 
be used in the description of the thermodynamical properties of associating 
fluids with arbitrary number of components, including the limiting case of 
poly disperse fluids. 
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1 Introduction 



Much of the progress achieved by the hquid-state integral-equation theories is due to 
the availabihty of the integral-equation approximations (lEA), which are amenable 
to the analytical solution for a number of the models of dense fluids and liquids. 
Since 1963, when Percus-Yevick approximation for the hard-sphere fluid was solved 
analytically [1, 2J, the analytical solutions were derived for a large variety of non- 
trivial Hamiltonian models (see [3], H] and references therein). During the last two 
decades substantial efforts have been focused on the development of the analyti- 
cally solvable lEA for the models of associating fluids [1]. Most of these studies 
were carried out in the frames of the product-reactant Ornstein-Zernike approach 
(PROZA) [SlElElElISlIinillllinilSlIIlIISlIiniinilig for the models, which 
combine hard-sphere interaction and sticky interaction due to a certain number of 
the sticky points located on the surface of each hard sphere. Different versions of 
the models, which describe dimerizing [3 El [9], polymerizing [TOl HIl [131 [El [IS] 
and network forming [TBI [T71 [TH] fluids, where investigated . More recently these 
studies were extended by adding a van der Waals attraction, modeled by a sum of 
the Yukawa terms [191 EQl [211 [22]. The properties of the corresponding models were 
studied using the analytical solution of the associating mean spherical approxima- 
tion (AMSA) [23]. In the particular case of the multi-component Yukawa dimerizing 
hard-sphere fluid with factorizable Yukawa coefficients the Yukawa part of the so- 
lution was reduced to the solution of only one nonlinear algebraic equation for the 
Blum's [211 [2S] scaling parameter F flUi [20]- In the limiting case of complete asso- 
ciation, when the system is represented by the multi-component mixture of Yukawa 
heteronuclear hard-sphere diatomics, solution of this equation represents full solu- 
tion of the AMSA. In this limit PROZA reduces to the 'proper' site-site theory [27] 
due to Chandler et al. [28] . 

In this article we extend solution of the AMSA obtained earlier [191 EQ] and de- 
rive explicit expressions for the thermodynamical properties of the multi-component 
Yukawa dimerizing hard-sphere fluid in terms of the F-parameter. 
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2 The model 



We consider M-component mixture of dimerizing Yukawa hard spheres of species 
i = 1, 2, . . . , M with diameters ai and densities pi. Each of the hard spheres has 
one sticky site placed on a surface. The pair potential of the model consists of the 
hard-sphere term, sticky site-site term and Yukawa term ^[J^ (r) , which was chosen 
to be of the following form: 

(3^(^\r) = -^e~'^ (1) 

where (3 = 1/kT, k is the Boltzmann constant, T is the absolute temperature and r 
is the distance between the centers of the spheres. 



3 Solution of the AMSA 

Solution of the AMSA for the model at hand was obtained earlier [19l [20] and we 
shall therefore omit the details here and present only the final expressions, which 
are needed in our derivation of the thermodynamics. 

AMSA consists of the two-density Ornstein-Zernike equation 

hij{k) = Cij{k) + '^piCii{k)cxihij{k), (2) 
I 

supplemented by the MSA-like closure conditions 

Cij{r) = E^e-^^ r > (Jij = {ai + cr,)/2 (3) 



h,j{r) = -E + -^5(r - (7^,), r < aij. (4) 

Znai 



Here hij{k), Cij{k), tij, a.i and E are the following matrices: 

hij{k) = \ ^ , Cij{k) =1 

,Ci^'jo{k) Ci^j^{k) j \Cjijo(^) '^hjii^) , 





tjjjj 



11 ai] 
















\ai j 




^0 
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where hi^j^{k) and Ci^j^{k) are Fourier transforms of the total hi^j^{r) and direct 
Ci^j^ir) correlation functions, respectively, ti^j^ = Tijgi^^j^Xtrtj), Qiojoi^^tj) is the con- 
tact value of the radial distribution function gi^jf^ (r) = hi^j^^ (r) + ^qo^/jo, Tij is the 
parameter which defines the strength of the sticky interaction and Ui is the fraction 
of non bonded particles. Here the lower indices a and (3 denote the bonding state 
of the corresponding particle and take the values (non bonded) and 1 (bonded). 
Fraction of non bonded particles ai together with association strength parameter 
Tij obey the mass action law (MAL) relation 

1 = f 1 + 2 ^ pj(Tijajti,j^ . (5) 



Solution of the AMSA was obtained using Baxter factorization method [29] with 
the general scheme of the solution based upon the version of factorization technique 
developed by H0ye and Blum [30l [31]. According to Baxter [29] the OZ equation 
(E]) can be factorized as 

S.,(|r|) = Q.,(r) -J2pi [ dr'Qu{r')c^iQji{r' - r), (6) 
I ^ 

J..(kl) = Q.,(r) + J2pJ ^"-'^y - ^l)«^Qo(^'), (7) 

I 

where T denotes the transpose matrix and the integrals Sjj(r) = 27i f^°° dr'r'cij{r') 
and Jij{r) = 2tx dr'r'hij{r') satisfy the following boundary conditions: 



Jjj(r) = Trr^E + Jjj, r < (Xij 

S,,(r) = E^e-^^ r > aij 



(8) 



Here Jij = Jij(O). From the analysis of the equations (Ej) and (JTj) we get 

Q^jir) = hj{r) + Uj] e{(Jij - r) + E^D,Je-"^ r > A^^, (9) 



where \ij = | (cTj — aj), Djj and E are the row vectors, i.e. Djj = {Di^j^, Di^^j-^) , 

E = (1, 0) . Vector Djj satisfies the following relation 

27r, 
z 



where Qij{k) = Sij {pjOtj) — 2tt J^ dr Clij{r)e^^^ and Kjj = {Kij,0). Expression 
for qij(r) in the interval Xji < r < (Jij is 

q,,(r) = ]-Wk, (r - a,,) (r - A,,) + E^^,- (r - a,,) + C,, (e^' - e-^-^ . (11) 



Here 



n \ n 



where 



S7jj(s) = 27rpjGij(s)Q;j, (^Gij(s) = dr'r'g. 
I I 



(12) 
(13) 

(14) 
(15) 

(16) 
(17) 



Here 'yij{z) satisfies the following set of the algebraic equations 



I 

and 



Tj = ^ piaiEcxitij, = ^ p,o-™, A = 1 - ^Ca, 



(19) 



1 — X — e 



>i a; 



1 - ^x - n + ix 1 e" 



One can see that all coefficients of the factor function Q(r) are determined by the 
set of unknowns Djj and 'y^jlz). These unknowns follow from the solution of the set 
of equations (fTUj) and ffT^ . 

Substantial simplification of the final algebraic equations representing the solu- 
tion of the AMSA occurs in the case of factorizable Yukawa coefficients, i.e. for 
Kij = Kdidj. According to Eq. (fTOj) now Djj can be written in the following form 



which gives 



/3i = ^^iE + Aia„ 



(20) 

(21) 
(22) 
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~ 27r / TT , \ - Ti . 27r _ 



Here 



vr 



(23) 

(24) 
(25) 
(26) 



where (^^{x) = (1 — e ^) /x. Next, making use of the symmetry property of the 
factor function, i.e. Qij{Xji) = Qjj(Aij), we have 



(27) 



where 



Xj = E^ [ciAi + die—'^'^J + a,Bj^o{zai)e-^"'\ (28) 
Equation ( l27ll enables us to introduce scahng parameter F via the following relation 

(29) 



a -2^X 

^3 — ^-J ' 



where D = J2k Pk^k^n-k'^k ■ Differentiating with respect to r and taking the 
limit r ^ we have 



D 



where 



R 



V 



(30) 



(31) 



Now all the unknowns of the problem can be expressed in terms of F, i.e. 

2A „ 
vr 



2A. 



vr 



A 



(32) 
(33) 



where nfJiT^ = -2A^j Y.kPk(^kthk^yk, {y = ?7,^, A), 



A,: 



1 + (po{zcXi)air' 1 + v9o(zcri)(TiF' 



2A 1 + ipo{zai)aiT' 
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^ 2^]^ (e^ - 1) - [2r + f C2 + ^ + 9^ 
' 6" (2r + f C2 + ^ + 2i]€) + A (z2 + ^Q^) (1 _ e«) ' ^ ^ 

f]^ = J]p, [azT,^ - 2// (1 - «^r,J] , 6^ = Y^p^ai {aiTf - yi) . (36) 

I I 

Finally, the nonlinear algebraic equation for F, which follows from (fTOj) . is 

(F)^ + zT + ttKD = 0. (37) 

Full solution of the problem requires solution of the set of equations formed by 
equations (ISl) and (1371) . The former equation needs as an input the contact values 
of the radial distribution function gi^jo- Corresponding expression follows from (I?!) 

Here we have used exponential approximation 



4 Thermodynamics 



Thermodynamic properties of the model at hand will be calculated via the energy 
route, which appears to be the most accurate for the MSA-type of the theories. Using 
standard expression for the excess internal energy in terms of the radial distribution 
functions, we have 

/3AE^ = -KJ2 PidiEcxiBj. (39) 

i 

Before proceeding to Helmholtz free energy calculations we will prove the following 
two useful relations 

dAE^ 

.+ . . \ 

V=const 
1 



_d (piaiPjajU^j^) _ 
dAE^ 



-^X X 



(40) 



dV 



Tip 



(F^ + zT) 



(41) 



pp=const 

where Pp denotes the set of all products piaiPjajti^j^. First of these relations can 
be derived by substituting Bj from ( l28l) into ( !39l) and differentiating it with respect 
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to piaipjajti^j^ with F held constant. We will skip these straightforward calcula- 
tions and proceed to the second of these relations. To prove relation (1411) we start 
eliminating X, from (EH]), (EO]) and ([29]) to get 



- EAi (^1 + ^zai + aiT^ - Errf,e-5-'^» - a^Er/ + ^ PfcXa^t^fc. (42) 

After substituting (|32l) into the above equation and performing some lengthy algebra 
we obtain the following matrix equation for Bj 



J2 BfcMfci 



1 TTCTj 

l + -^a. + a.r + — G^E-r..^^^, 



2n ^ f 1 \ , 



^Pfcttfe^njiAfc (43) 



^ fc ^ ^ / k I 



where J = (0, 1), M^, = e^-""^ (1 + {zau) crkT) Pfc„ 

(Ti {l + \zai + OiV + ff-C2} V^o {zfJi 



Xi 



1 + V9o [zai] aiV 



(44) 



and P is the Jacobi type of the matrix, i.e. P^i = ^ifcl + c^^dj + e^fi. Four vectors 
c, d, e and f , that form the Jacobi matrix, are 



Cj = Piii (1, ttj) , di = (Ti, 

















^0 





Corresponding equation for 
respect to V 



9B, 



follows from (H3ll upon its differentiation with 



gBfc 

9r 



M 



ki 



-Xj — J ( — y^Pfcctfctiifci— -^fco 



(45) 



Taking derivative of the both sides of equation (!39|) with respect to F, using expres- 
obtained from the solution of the set of equations and taking 



sion for 



9B, 
dV 



Pp 



into account equation (!37|) . we recover relation (14T!) . Inverse matrix M^^^ together 
with r derivatives, which are used in our calculations, are given in the Appendix. 

We start our derivation of the expression for Helmholtz free energy with the 
following standard thermodynamic relation: 

^ (/?AA) = f3AE. (46) 

Integrating this equality by parts and using the fact that all thermodynamic quan- 
tities according AMSA depend only on one parameter F, we have 

PAA"- = PAE"- - £ dT'p'^, (47) 

where AA^ represent Yukawa contribution to Helmholtz free energy. We start with 
the system of dimerizing hard spheres and charge it by the Yukawa charge di up to 
the current conditions. Full derivative under the integral in f jTfl) can be expressed 
in terms of the partial derivatives giving 



Pp ij 



dAE 



Y 



dpiaiPjajt 



(4J 

Integrating the second integral in (HHI) by parts and using (p7|) and (HTI) . we get 

1 /F^ F2~ 
(3AA^ = (3AE^+'-('- + zL 

TT \ 6 I 

+ KY^PiaiPjajti^j.Xi^Xj,^-^ j dV piaiPjajU.j, ^° • (49) 

ij ij 

Since dtj-^ki/dT = tj^^kidln^gj^^ko) /dT and due to the exponential approximation 
flH^ and MAL relation it is straightforward to show that 
d PA^'^^ ^ d KX,,Xj, 

where 

p/^jlMAL = J2\nai + J2 PiCiipiajU.j, . (51) 

i ij 

The final expression for Helmholtz free energy in excess to Helmholtz free energy of 



dimerizing hard-spheres system is obtained combining (1491) and (jSOll 

13 AA^ = I3AE''+^ ("y + z^-^\ +K A«iP,«,tn.l^^o^.o+/?A^*'^''-/5A<^^ 

(52) 



where 



(53) 



In the case of the reference system represented by the multicomponent hard-sphere 
mixture we have 



1 /r^ 

PAA = /3AE + -l— + z—] +Kj2Pi(^iPj(^jthh^A + PAA^^"^, 



(54) 



with AA being Helmhohz free energy in excess to the hard-sphere Helmholtz free 
energy. Corresponding expression for the excess entropy AS is found differentiating 
( IMl) with respect to the temperature 



k 



AS = -— ( V + -^V ) - kB/3AA 



[MAL 



ksK ^ piaiPjajti^j^Xi^^Xj^y (55) 



Similar as in the earher studies [32| [33] the scahng parameter F of our theory 
minimizes the excess Helmholtz free energy 



4aa = o. 



(56) 



Differentiating with respect to the density of one of the components pz, we get 
expression for the chemical potential 



dAE 



dpi 



+ - (r^ + zV) 

TT ' 



dpi 



dA 



MAL 



dpi 



d piOf-i Pj oijti-^j-^ 
dpi 



-^io-^jo ~^ Pi^iPj^j^iiji 



dXi^Xj^ 




dpi 





(57) 



Using (I5l) and (l38i) . it can be shown that 

P = liaai- 2_jPi^iP3^i^hh y^i. 

According to ( 140|) and ( 14T1) we have 



^ln4^(a.,) 
dpi 



+ K 



d {Xi^Xj^ 



JO J 



dAE 

dpi 



--{T + z) 



IT 



dr_ 

dpi. 



-K^Xi^Xj^ 



dpiaipjajU^j, 
dpi 



-1/3 



dpi 



dAE 
. dpi 



1/3 ^J 

The latter two expressions, when substituted into (jSZ]), yield the following simple 
expression for the chemical potential 



PAfii = (3 



dAE 
dpi 



Pp,T,0 



+ In - ^ PiaiPjUjU^j^ 



1-3 
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a In 4^ (a,,) 



dpi 



(58) 



Expression for 



dAE 

dpi 



was obtained using (l39|l and (l28|) . it reads 



/3 



dAE 



-Kd, 



-zai/2 



-K ^ pidi 



-Z(Ti/2 



- Jpp,r,/3 ""'ai(po{zai) 
9Ai 



X. 



cO 



+ KAidi 



g-Z(7;/2 



+ Kd1 



-zai 



where derivatives 



dAi 
dpi 



aiifo (zai) 



dpi 



+ 



dpi 



TT V 



9?7 



dpi 

Z(Ti/2 



(59) 



and 



are given in the Appendix. Finally for 



the excess pressure AP one can use the following standard relation: 

/3AP = /?5^p,A/i,-/5AA 

i 

Remarkable fact is that chemical potential and pressure are independent of 
and thus we don't need to solve any equations to obtain this derivative. 



(60) 



ar 



dpi 



5 Summary and concluding remarks 

In this paper we consider multi-component dimerizing Yukawa hard-sphere fluid. We 
present explicit analytical expressions for Helmholtz free energy, chemical potential, 
entropy and pressure of the system in terms of the Blum's scaling parameter F, 
which follows from the solution of the AMSA for the model with factorized Yukawa 
coefficients. In the latter case solution of the AMSA reduces to the solution of only 
one nonlinear algebraic equation for F. This feature enables the theory to be used 
in the description of the structure and thermodjTiamics of associating fluids with 
arbitrary number of components, including the limiting case of polydisperse fluids. 
We are currently studying the effects of polydispersity on the phase behavior of 
the polymer fluid combining the theory proposed here and dimer thermodynamic 
perturbation theory for polymers [311 135] . 



6 Appendix 

We present here expressions for the elements of the inverse matrix M^^ and deriva- 
tives, which are needed to prove relation (HTl) and appear in the expression for the 
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chemical potential (158|1 



'dvi' 




2A 


9^ . 


Pp 








2A 


dV 


pp 





where y takes the values ^,77 or A 

' dxi 

or 



2A ^ v-^ 2A ^ 

fc 



2A 2A 

"S«Xi ~l Si- 



TTCT,; 



(61) 
(62) 



(63) 



_ ZCTj 

e 2 



1 + V^o (zai) aiT 
Az^ (1 - e«) 29 



^ ^e^fi + ^c^t; + 



nS 

where denominator 



fcdi 



^ = e" (2r + ^C2 + z + 2Q^) + 7 + J^") (1 - 0^) • 



(64) 



(65) 



, dpi 

dAi 
dpi 



S 

(2^]5 + 2r + ^ + ^) (^^7]^ + x,o) - (1 - e«) Li 



where 



= 2r^V(l 
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+ 



2A 



6A y ' 3A + + +2X,o. 



(66) 
(67) 



(68) 
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